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The "s-f model," also known as the ferromagnetic Kondo 
lattice, contains a description of band electrons coupled to lo- 
calized spins which is an appropriate description of the mag- 
netic part of the low-energy physics of Gd metal. Here the 
model is used to estimate the lifetime broadening of the mi- 
nority spin component of the surface electron band in ferro- 
magnetic gadolinium metal at temperatures below the Curie 
temperature. The low temperature result 1/r «0.1 eV agrees 
nicely with a measurement by Fedorov et al. 

I. INTRODUCTION 

Fedorov et al. W have recently measured the line-shape 
A{k, uj) of photoemitted electrons in the ferromagnetic 
metal Gd. For photoelectrons associated with a pho- 
toholc in the (001) surface band, they find somewhat 
different life-time broadening depending on whether the 
emitted electron has up (majority) or down (minority) 
spin. They interpret the source of lifetime broadening 
to be electron-phonon scattering for the majority spin 
component of the photohole, and electron-magnon scat- 
tering for the minority spin component of the photohole. 
Although the arguments given by Fedorov et al. seem 
perfectly sensible, nevertheless, this interesting diversity 
suggests a need for theoretical inquiry. The magnitude of 
electron-phonon scattering in Gd has previously been 
estimated, with results roughly agreeing with the as- 
signed majority spin equilibration rate. Electron-magnon 
scattering has not previously been estimated. 

Here I argue that the "s-f" or "ferromagnetic Kondo 
lattice" model allows reasonable estimates without free 
parameters. I make an extreme model for the nature 
of surface electron and surface magnon states: surface 
electron states have amplitude 1 on the top layer and 
zero elsewhere, while magnon states at the surface || are 
simply the bulk Bloch states, ignoring surface boundary 
conditions. Using this model, and the measured mass 
to* w 1.2 of the surface hole band, the zero tempera- 
ture equilibration rate of minority spin holes more than 
25meV from the top of the hole band is predicted to be 
1/r =0.10 eV, agreeing with the experiment. 

II. FORMULA FOR RELAXATION RATE 

The generic Hamiltonian for the "s-f" or "ferromag- 
netic Kondo lattice" model couples electron bands kn 



with energy e(kn) (independent of spin, so far) to local- 
ized spins Si located on atoms at lattice sites Rf. 

TL = e(kn)c* (kruj)c(kno) 

knu 

— J Si ■ C 1 (ima)<j a pc(im/3). (1) 

ima/3 

The electron bands derive from the outer atomic or- 
bitals \ima > with wavcfunctions ip m (r — Ri)Xa where 
Xa is the spin part. This model should describe the low- 
energy spin-related physics of a metal like Gd, with three 
conduction electrons per atom in orbitals derived from 
atomic s and d states of Gd. These three orbitals have 
a magnetic interaction with the well-localized S = 7/2 
half-filled 4/ shell. The exchange parameter J comes 
from the atomic Hund's rule trying to keep "electron" 
spins parallel to "core" spins. The same Hamiltonian, 
in the J — > oo limit, is known as the "double exchange 
Hamiltonian" Q and is very popular right now || for 
discussions of hole-doped LaMnOs. 

The zero temperature phase diagram of the s-f model 
(for a single s band) as a function of filling and J/band- 
width ratio has been computed approximately || . Ferro- 
magnetic order occurs over a wide range of parameters, 
with a Curie temperature (T C =292K for Gd) proportional 
to J. Electron bands acquire a spin splitting proportional 
to J. Lindgard et al. JjJ used the model to calculate 
(in random-phase approximation) the spin-wave disper- 
sion, which was measured for Gd by Koehler et al. ||, 
and which has been studied using spin-density functional 
theory by Perlov et al. ||. 

In lowest-order spin-wave approximation, we replace 
the Fourier-transformed spin operator Sq in Eq. ([!]) by 
spin-wave creation and destruction operators Oq and oq 
using Sq z = S6q.q - a} Q a Q , (S Qx + iS Qy )/2 = S^ = 

V2Sciq, and {Sq x -iSq v )/2 = Sq = V2Sa Q . The Q = 
term in lowest order gives spin splitting 2 JS, lowering the 
energy of bands with spin parallel to the localized spin 
Si ~ Sz and raising the other bands equally. The S^ 
terms give rise to spin-flip scattering events. 

In this note I estimate the rate 1/r at which a single 
out-of-equilibrium hole in a surface band relaxes back to- 
ward equilibrium by spin-flip processes. The rate can be 
found from H/t — — 2ImE, where the leading Feynman 
diagram for the self-energy £ is shown in Fig. |l|. Only 
the one-magnon process is considered. This approxima- 
tion can be questioned on the grounds of fallibility of the 
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"Migdal approximation" for electron-magnon processes 
|^o| . However, it is a proper first estimate, the only one 
which can be reliably computed, and seems to me un- 
likely to make a large error when T « T c . Two magnon 
processes have been considered by Lutovinov and Reizer 
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FIG. 1. Feynman graph for electron self-energy from elec- 
tron-magnon coupling. The double solid line is the renormal- 
ized electron Green's function; the wavy line is the renormal- 
ized magnon Green's function; circles are the effective elec- 
tron-magnon interaction matrix element. 

The actual states (fc, n) of interest are surface states of 
spin primarily but not exclusively up. The probability 
p(k,n, |) (or p(|) for short) that the spin state is up, is 
less than one. The amplitude of the up-spin component 
of the wavefunction is y/p("\) and is close to 1. The cor- 
responding probability 1 — p(^) = p(k, n, \) (or p([) for 
short) that the spin state is down, is non-zero for two 
reasons. First, the spin-orbit interaction is not small and 
mixes spin states. Second, /-spins on Gd atoms may 
deviate from perfect alignment by quantum and ther- 
mal fluctuations, and the conduction states are locally 
locked by Hund's rule in the same spin orientation as the 
/-spins. The second process is a renormalization of elec- 
tron bands by magnon processes. A combination of the 
two effects is seen experimentally as a small minority 
spin component in the photoemitted electron. This com- 
ponent is sometimes referred to as a "shadow band" and 
has received recent theoretical treatments |l2| . 

An elementary derivation of the "Migdal" result fol- 
lows from the standard "Golden rule" rate equations 
of Boltzmann theory. Suppose the down-spin compo- 
nent of state kn has population p(l)F(k, n, }). If this 
deviates from the equilibrium Fcrmi-Dirac population 
p(l)f(k, n, |), then it will evolve back toward equilibrium 
according to 



dp(i)F(k,n,j) 
dt 



2n 



\M b (kn i,k + Qn' T)| 2 



Qn'b 



x {L(emission) + L( absorption)} . (2) 

L (emission) = S(e — e' — u>') 
x [(N' + 1)F(1 - F') - N'F'(1 - F)\ (3) 



Here M b () = J[(2S/2N)p(i)p'(f)} 1 / 2 is the matrix ele- 
ment for the process (fc, n, }) scattering to (fc + Qn' T) 
by emission of the magnon (—Qb) or absorption of the 
magnon (Qb). The factor l/\/2N, where N is the number 
of unit cells and 2N the number of atoms in the sample, 
comes from normalization of the spin-wave eigenvector. 
There are two branches of spin waves, with amplitude 
±I/\/2 on each atom. This is discussed in the Appendix. 
A shorthand is used that e and p(l)F stand for the en- 
ergy and occupancy of the quasiparticle state (kn |), e' 
and p(1)F' stand for the energy and occupancy of the 
quasiparticle state (fc + Qn' |), lo and N stand for the 
energy and occupancy of the magnon state (Qb), and 
lo' and N' stand for the energy and occupancy of the 
magnon state (—Qb). At temperature T > 0, depletion 
of an excess population toward equilibrium occurs both 
by emission and absorption of thermal magnons. Each 
process (emission or absorption) has a time-reversed pro- 
cess which enhances the population, the "scattering in" 
terms with opposite sign. In thermal equilibrium, scat- 
tering out and in occur at equal rates. This "princi- 
ple of detailed balance" guarantees that the two parts of 
L (emission) cancel each other when the distributions N 
and F become the equilibrium distributions n and /, and 
similarly for L (absorption). 

Now make the assumption that all quasiparticles are in 
equilibrium except for a particular state (kn I) of inter- 
est, whose population (F) deviates from equilibrium (/) 
by 5F(kn J,). Then the rate equation || takes the form 



dF(k,n,l) 



dt 



= -5F(kn l)/T(k,n, [) 



(5) 



UN ^ 

Qn'b 

x {5(e - e' - u)[n+ 1 - /'] + 6(e - e' + u)[n + /']} (6) 

Except for the factor £>'(}) = p(k + Qn' })> this magnon- 
limited scattering rate is a perfect analog of the usual 
phonon-limited quasiparticle relaxation rate from Migdal 
theory. Eq. @ can equally well be derived by evaluation 
of the Feynman diagram Fig. [I], analytic continuation 
of Matsubara frequencies ilo v to complex frequency z, 
and use of 1/r = — 2ImS(z — > e + i8). The energies 
e and ui used to evaluate Eq. (||) are the experimental 
quasiparticle and magnon energies, which is mandated 
by the occurrence of the full, not bare, Green's functions 
in the self-energy graph of Fig. [|. 



III. EVALUATION OF RELAXATION RATE 



L(absorption) = 8(e — e' + lo) 
x [NF(l-F')-(N+l)F'(l-F)} (4) 



Reasonable simplifying assumptions yield a "zero pa- 
rameter" estimate of Eq. (||) for 1/t. First, we specialize 
to a state (fc, n) which is in an occupied surface band. 
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This state at low T is dominantly spin up, but as seen ex- 
perimentally, it has a spin down component, p{[) w 0.13. 
The Green's function G and self-energy are 2x2 matri- 
ces in spin space. Eigenstates of the 2x2 matrix G _1 
are the spin-split quasiparticles, of which only the lower 
(primarily majority spin) state is relevant. The energy of 
this state will be denoted e k , with no band or spin index 
needed, and only a two-dimensional wavevector k rather 
than a three-dimensional k. Spin-resolved photoemission 
selectively depopulates a single spin component of this 
state. When the photoemitted electron has spin down, 
there is a hole in the down-spin component of the of the 
majority spin hole band. This component (with distri- 
bution function F{k [)) decays to equilibrium with rate 
1/t(/c,|). We assume that p(k J.) = p(J.) is independent 
of k, and similarly p(k + Q f) = p'(t) = 1 - p(i)- Then 
we have 



h/r(k,i) = ^(2JSfV 



(7) 



where V is the relevant density of decay channels. At 
T = 0, V is 



V 



-Yn 



y/3a 2 c/2 
(2tt) 3 



x / d 2 QdQ z {S(e k - e k+Q - oj(Q, Q z ))6(+e k+Q ) 

+S(e k - e k+Q + u){Q, Q z ))6(~e k+Q )} (8) 

We assume that the spin waves retain their bulk charac- 
ter. Therefore, although electron energies depend only 
on the two-component k, the spin-wave energy u(Q, Q z ) 
depends on all three components of wavevector, the z 
direction being normal to the surface. The first delta 
function in Eq. (|J) is a process where an empty state 
k + Q lying above the Fermi surface is filled from a state 
k which scatters into it by spin-wave emission. In other 
words, it refers to decay of an electron lying in a state k 
which is above the Fermi energy by at least a spin-wave 
energy. The second delta function is a process where the 
filled state k + Q lying below the Fermi surface scatters 
into an empty state k by spin-wave emission; in other 
words, it describes decay of a hole in the state k (below 
the Fermi energy by at least a spin- wave energy.) It is 
only this second process which is seen in the photoemis- 
sion data. Let us also assume that the photohole state 
k lies below the Fermi energy by at least 25 meV, the 
maximum spin-wave energy, in which case 0(— efc+g) is 
guaranteed to be 1. The Q z integration then gives 



2tt 



7r/c 



t/c 



dQJ(e k - e k+Q + uj{Q, Q z )) = v{k, k + Q)/s(Q) 



(9) 



where s(Q) is the normalized spin-wave slope 



s(Q) 



2n 

c 



dQ z 



(10) 



and v(Q) is the number of spin-wave states (Q, Q z ) with 
fixed Q, but any value of Q z , which conserve energy, i.e. 
which satisfy oj(Q, Q z ) = uj* = £fe+Q — efe- Examining the 
measured dispersion curves || , the normalized slope can 
lie between and a maximum which is not very different 
from a; max =25 meV. If efe+Q — e k is not greater than w max , 
then there is a fairly good chance (something like 50% 
probability) that there are 2 values, ±|Q Z |, such that the 
spin-wave state (Q,Q Z ) obeys energy conservation. In 
other words, v{Q) can be expected to take the value 2 or 
with about equal probability provided < e k +Q — e k < 
w mal , and is definitely elsewhere. 
The remaining integral is 



V = 



y / 3a 2 /2 
(2nf 



d 2 Qv(k,k + Q)/s(Q). 



(11) 



This integral should be evaluated using the correct dis- 
persion relations for the surface electron state and for 
the spin- wave states, with allowance for boundary condi- 
tions and wave-function amplitudes altering the matrix 
elements. However, this requires a large "first-principles" 
calculation of uncertain reliability. Therefore, a slightly 
cavalier estimate is in order. Assume that the surface 
state has parabolic dispersion e k = — eo — h 2 k 2 /2m* . Ex- 
perimentally this state is seen jjj to disperse downwards 
in energy with effective mass m* ~ 1.2m, where m is the 
electron mass. 

The inequality < t k +Q — e k < w max is obeyed in an 
annular region shown in Fig. 0. The area of this region 
is n(k 2 -k' 2 ) and (h 2 /2m*)(k *- k' 2 ) = w max . Inside this 
region we assume v has an average value of 1 and the 
normalized slope has an average value w max /2. Then the 
density of decay channels is 



V 



V3a 2 /2 2m* 
0.57(eV) _1 



%/3 



ma 



2n h 2 



(12) 



This is e-independent because of the 2-dimensionality of 
the surface band. A more careful treatment, yielding 
exactly the same result, is in the Appendix. 

Finally we evaluate the decay rate Eq. (Q) by choosing 
2JS to be the spin splitting of the surface state, ~ 0.65 
eV |fL3l , S to be 7/2, and the fractional up-spin probabil- 
ity p' ~ 0.87 §]. These assumptions yield %/t - 0.10 eV 
for all photohole states k which are not too much closer 
than 25 meV to the top of the surface state band. For 
states closer to the top, the decay rate should diminish 
because of reduction of the number of decay channels V. 
It is also assumed that T is fairly low. At higher T, the 
rate is enhanced by thermal spin-wave population, but 
diminished by diminishing fractional up-spin probability 
p'(T). The estimate %/t ~ 0.10 eV agrees with experi- 
ment 0] to greater precision than the uncertainties of the 
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model. The calculation supports the interpretation that 
the source of minority spin line-broadening is spin-flip 
decay. 




FIG. 2. For a given hole state k, the annular region shows 
possible lower energy holes (higher energy filled states) which 
the hole can scatter into by spin-wave emission. The outer 
circle denotes states with energy — eo — h 2 k 2 /2m* , and the 
inner circle denotes states higher in energy by w ma x. 
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the second atom in the two-atom basis of the hep crystal 
structure. The frequencies and density of states calcu- 
lated from this model fit experiment [0 well, and are 
shown in Fig. 0. 

Notice that the two branches u>± are degenerate at the 
top of the Brillouin zone Q z = tt/c (directions A — ► H — > 
L —>■ A in Fig [|) These branches can be "unfolded" 
into a zone twice as large in the z-direction; u>+ is the 
extension of U- under the mapping LU + (Q Xl Q yi Q z ) = 
u)-(Q x ,Q y ,2ir/c- Q z ). 



> 25 

CD 
O 

§ 15 

cr 


I 10 

> 

5 5 
I 



K M 



r a 



H 



A D(co) 



FIG. 3. Spin wave dispersion in Gd, calculated from 
Lindgard's parameters [j?]. The magnon density of states 
is shown in part (b) of the panel. The thin line represents 
the dispersion curve obtained by truncating farther than first 
neighbors in the z direction. 



APPENDIX: 

Here I try to illuminate some of the approximations in 
the model and also to demonstrate the surprising relia- 
bility of the estimate 



VV(T)to* (2JSa 



(Al) 



which comes from combining Eqs. (Q) and (|l^). Spin 
waves in Gd have been fitted Jl4| with the model 



Vii(Q) V 12 (Q) 
V 12 (Q) Vu(Q) 



1/V2 \ 



■■v±(Q) 



l/y/2 \ 
±1/V2j 

(A2) 



where the spin wave eigenfrequencies are uj±(Q) = 
Vn (Q) ± V\2 (Q)- The elements of the matrix are defined 
by 

V n {Q) = SVii(O) - Ju(Q) + Ji 2 (0)] 

V 12 (Q) = -SJ 12 (Q) (A3) 

J n (Q) =ReJ] J(R)e^ A 

R 

J 12 (Q)=ReJ2j(R + ?y^ +r) (A4) 

R 

where J(R) are fitted exchange coupling constants, R 
runs over the translation vectors of the hexagonal lat- 
tice, and f with z-component c/2 gives the position of 



An accurate analytic approximation for Eq. (|ll| ) 
i s p ossible provided Lindgard's fitting parameters, Eq. 
( A4 ) are slightly simplified by truncating off farther than 
first neighbor planes in the z direction. The alteration 
of uj(Q) caused by this approximation is shown in Fig. 
||. The change is fairly small, and could be largely 
compensated by a further tuning of the nearer neighbor 
couplings. With this approximation, the one extended 
magnon branch uj = u>- has frequency given by 



u>(Q) = a(Q x , Q y ) + b(Q x , Q y ) cos(Q z c/2). 
Then the normalized slope s(Q) is given by 



s(Qx,Q v ) = Tr\/b 2 - (u> - a) 2 



(A5) 



(A6) 



The remaining integral Eq. ([11]) can be considered an 
integral over cPk' running over the annulus of Fig. 0, 



V 



V3a 2 / 2 



k 2 

d(j)k dk — - — 
fci n^b 2 



(u - a) 2 



(A7) 



where uj = e(k) — e(fc') is the energy of the spin wave 
with wavevector (k x — k' x ,k y — k' y , Q z ) which scatters the 
hole out of state k into state kl . For a given two- vector 
(Qx, Qy) — k — kl , one searches over Q z to find whether 
there is an energy conserving solution. Either there are 
no energy conserving spin-wave states, or else there are 
v = 2 such solutions at ±Q Z . Consider a path in k'- 
space at fixed azimuth <f> shown as a dotted line in Fig. 
|2|. The outer circle is states k' which are degenerate in 
energy with the starting state k. Unless k' is the same 
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as k, there is no zero-energy spin wave which can couple 
these states. Moving down the dotted line to higher en- 
ergy electron states which can fall into the hole at k by 
magnon emission, one finds the state labelled 1 which has 
just enough energy difference that an energy conserving 
magnon transition is found. The magnon has Q z = 
because it is the least energy magnon allowed to cou- 
ple on the dotted line. Moving farther down the dotted 
line, one comes eventually to the state labelled 2 which is 
the highest energy electron state which can fall into the 
hole in state k by magnon emission. The magnon has 
Q z = 2-k/c because it is the highest energy magnon. For 
both these extreme states, u) — a ± b, the slope s is zero 
(the magnon energy is quadr atic in Q z near Q z = and 
2ir/c.) The integrand of Eq. A7 diverges at the two end 
points. However, it is an integrable divergence, and in 
fact, almost exactly independent of the variables a and 
b. We observe that in going from states 1 to 2 along the 
dotted line, the components Q x and Q y are not changing 
much, which allows us to set a(Q x , Q y ) and b(Q x , Q y ) to 
constants during the fc'-integration. Then wc have 



V = 



V3a 2 /2 rrf 
(2tt) T 2 

\/3a 2 m* 



-a+b 



de'- 



2/tt 



- a -b y/b 2 - (e - e' - a) 2 



(A8) 



This is exactly Eq. (jT|). 

Finally it is appropriate to mention the hidden as- 
sumptions about surface and bulk states. It is implic- 
itly assumed that the electron surface state has a wave- 
function of unit amplitude on the surface layer and zero 
amplitude elsewhere. This kind of state will be ab- 
solutely insensitive to the z-component of the magnon 
wave- vector, and the electron-magnon matrix element 
will be J[(2S/2JV>fl>(i)] 1/2 for all k and k' . Sup- 
pose instead a surface electron wavefunction of ampli- 
tude l/-\/2 on each of the top two layers. This state 
will couple to the lj- branch of magnons with the full 
matrix element J(2S/2N) 1 / 2 , but will not see the 
branch. Conversely, a surface wavefunction which has 
amplitude ±l/v2 on the top two layers, with a sign 
change, will couple only to the w+ branch and not to 
the lj- branch. It is not reasonable that the net cou- 
pling to magnons should depend much on the depth or 
details of the surface electron wavefunction; therefore we 
should ask where has the missing magnon coupling gone 
when the electron wavefunction extends two layers down 
instead of one. The answer is in the orthogonal electron 
wavefunction on the top two layers, with opposite phase 
relation between layer 1 and 2. Since by assumption there 
is only one surface state, the orthogonal state is not an 
eigenstate but a superposition of bulk states. The miss- 
ing magnon coupling is from the surface state into these 
bulk states. For the actual Gd surface state, how much of 
the magnon-induced scattering is to bulk and how much 
is to surface states is an unknown element. The extreme 



model used here hides this problem. The justification is 
belief that the net scattering probability should have a 
tendency to be conserved, i. e. to be weakly dependent 
on depth. Similarly, we have not asked what is the na- 
ture of the magnon states near the surface, but instead 
assumed that we can use bulk magnon states. Instead, it 
might be that a surface band of magnon states grabs all 
the spectral weight. Then details would be quite differ- 
ent, but over-all coupling strength should be similar. 
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